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FUSE OBSERVATIONS OF THE HD MOLECULE TOWARD HD 73882 
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ABSTRACT 

The Lyman and Werner band systems of deuterated molecular hydrogen (HD) occur in the far UV 
range below 1200 A. The high sensitivity of the FUSE mission can give access, at moderate resolution, to 
hot stars shining through translucent clouds, in the hope of observing molecular cores in which deuterium 
is essentially in the form of HD. Thus, the measurement of the HD/H 2 ratio may become a new powerful 
tool to evaluate the deuterium abundance, D/H, in the interstellar medium. We report here on the 
detection of HD toward the high extinction star HD 73882 (£b-v = 0.72). A preliminary analysis is 
presented. 

Subject headings: ISM: abundances — ISM: clouds — ISM: lines and bands — ISM: molecules — 
ultraviolet: ISM 



1. INTRODUCTION. 

It has long been recognized that the primordial abun- 
dance of deuterium represents the most sensitive probe of 
the baryonic density Qb of the Universe (see, e.g. Schramm 
& Turner 1998; Olive, Steigman & Walker 1999). More- 
over, the abundance of deuterium at any epoch is a lower 
limit to its primordial abundance since deuterium is solely 
destroyed in stars of any mass. Therefore, D/H is also an 
efficient tracer of the universal star formation rate. How- 
ever, the evolution of the deuterium abundance from zero 
to solar metallicity is still unclear. 

Measurements of the atomic (D/H) ratio have been per- 
formed in different astrophysical sites, namely in moderate 
to high redshift quasar absorbers, in the presolar nebula 
and in the local interstellar medium (for reviews, see e.g. 
Ferlet & Lemoine 1996; Linsky 1998; Vidal-Madjar et al. 
1998a; Lemoine et al. 1999). These studies indicate that 
D/H may vary within the local interstellar medium by a 
factor as high as ~ 2 over spatial scales of few tens of par- 
sees (Vidal-Madjar et al. 1998b; Jenkins et al. 1999; Son- 
neborn et al. 2000), while presolar and quasar absorbers 
D /H abundances are limited by the existing scatter in the 
results. 

Deuterated molecules are another means of estimat- 
ing the deuterium abundance. To date, over 20 single 
D-bearing species and two doubly deuterated molecules, 
D 2 CO and ND 2 H, have been observed at radio frequen- 
cies both in cold interstellar dark clouds and in warmer 
star forming regions (see e.g. Roueff et al. 2000). How- 
ever, chemical fractionation takes place in cold regions and 
mantle desorption of grains are often invoked in star form- 
ing regions. Deriving accurate deuterium fractional abun- 
dances using these molecules is therefore very difficult. 



Recently, the R(2) transition at 37.7 /im of HD has been 
detected with ISO in giant planets (Feuchtgruber et al. 
1999) and the pure rotational J=1^0 line at 112 /im to- 
ward the Orion Bar (Wright et al. 1999). Bertoldi et 
al. (1999) have also detected the excited rotational line at 
19.43 /im of HD J=6^5 in Orion KL, a molecular outflow 
region. Although ISO thus opened the sky to HD emis- 
sion, the derived column densities depend strongly on the 
modelling of HD excitation and on extinction corrections. 
With Copernicus, H 2 and HD molecules were observed in 
absorption in the ultraviolet in diffuse interstellar clouds 
such as that toward C Oph (Wright & Morton 1979). How- 
ever, the low HD/H 2 value found (fewx 10~ 7 to fewx 10~ 6 ) 
reflects the mostly atomic nature of these diffuse clouds; 
to determine the D/H ratio from these data requires a de- 
tailed model of the formation and of the destruction of 
HD. 

Because of its high throughput, the Far Ultraviolet Spec- 
troscopic Explorer (FUSE; Moos et al. 2000) can give ac- 
cess to denser molecular clouds within which one might 
expect H and D to be essentially in their molecular form. 
If so, no chemical fractionation correction will be required 
and accurate deuterium abundances could result directly 
from the measurement of the HD/H 2 ratio. In Section 2, 
we present FUSE observations of HD in the translucent 
cloud in front of the star HD 73882, along with the data 
reduction and analysis. A companion paper by Snow et 
al. (2000) describes the H2 observations. A preliminary 
discussion is given in Section 3, while some conclusions are 
drawn in Section 4. 

2. FUSE OBSERVATIONS OF HD 73882. 



Twenty-seven years after Copernicus, FUSE was suc- 
cessfully launched on June 24, 1999 from Cape Canaveral. 
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It has a sensitivity of about 10 times that of Copernicus, 
in the wavelength range from slightly below the Lyman 
limit (905 A) to 1187 A (Moos et al. 2000). The higher 
throughput comes about because of the use of a guidcr 
that can be used for fainter stars and of the use of an ar- 
ray detector instead of a scanning photomultiplicr tube. 
Many absorption lines from different rotational states of 
H 2 and HD are available within that range. 

The bright early type (08.5V) star HD 73882 at about 
925 pc has been the first target shining through a translu- 
cent cloud, with a reddening i?B-v=0.72 or A v =2.44 (note 
that £^b-v=0.38 and A v ~1 were the highest values ob- 
servable by Copernicus) , to be observed with FUSE in the 
time-tagged mode, for a total integration time of ~ 7 hours 
between 30 and 31 October 1999. At this time, the spec- 
trograph was still unfocused, resulting in a spectral resolv- 
ing power R=A/AA ~ 14000. Unfortunately after a short 
time the two short-wavelength SiC channels became mis- 
aligned and failed to accumulate on the target. However, 
for our very reddened target, most of the information con- 
cerning HD is expected to lie above 1000 A, in the LiF 
channels. 

The spectra were extracted from the 11 sub-exposures 
with the version 1.6 of the pipeline, and the electronically 
induced detector bursts were subtracted by screening of 
the events with IDL tools. The processing includes data 
screening, thermal drift correction, geometric distortion 
correction, Doppler correction to heliocentric wavelengths, 
dead time correction (data loss due to overflow of photon 
events) and wavelength calibration, but not correction for 
astigmatism. In order to obtain the best resolution, the 
calibrated spectra were then manually coadded over small 
wavelength windows, using both exposure time and error 
weighting. The resulting signal to noise ratio is ~ 20 per 
resolution element in the continuum near the strongest 
available HD lines. 

Fig. 1 shows the final recorded spectrum. Broad 
H2 bands are clearly seen (see also Snow et al. 2000), 
together with some CO bands and several other narrow 
interstellar features including atomic species such as Ci, 
Ol, Nl, Fell and Pn. One can note that the heavily sat- 
urated H2 lines reach a zero flux, indicating an extremely 
low level of scattered light in the FUSE instrument. 

We have performed a preliminary analysis of these data 
by using a simultaneous profile fitting of all detected lines 
from any molecular and atomic species included in a recent 
updated compilation by Morton (private communication). 
In a first step, a single absorbing component was assumed 
(and found at a coherent radial velocity for all lines). It 
has to be noted that most of the lines arc saturated so that 
the column densities are highly sensitive to the value of the 
intrinsic broadening parameter - b- value - which combines 
thermal and turbulent widening. Furthermore, the uncer- 
tainties in the instrumental line-spread function make it 
quite plausible that large amounts of cold absorbing gas 
might be " hidden" . Nevertheless, the dominant J = and 
J = 1 levels of the H 2 lines lie on the damped part of the 
curve of growth and the total H 2 column density is al- 
most independent of the assumed 6-value and the velocity 
structure of the line of sight. We find log 10 -/V(H 2 )t t = 
21.2 ±0.2, in agreement with Snow et al. (2000). 

Absorption from R(0) lines of HD are also detected to- 



ward HD 73882, four examples of which are presented 
in Fig. 2, arbitrarily set to km s" 1 radial velocity. 
This is a ~ 5<r detection, as shown in the upper panel 
of Fig. 3, in which is plotted the A% 2 of the simultane- 
ous fits of all detected HD lines (except the one blended 
with Ovi near 1031.9 A) as a function of the HD column 
density: when 7V(HD) decreases, A% 2 tends to nearly 25. 
The HD lines, as well as the CO and the high J level 
H 2 ones, lie on the flat part of the curve of growth where 
AT(HD) heavily depends upon the assumed effective b- value 
(note that this makes useless the use of equivalent widths) . 
In that case, the assumption of a single absorbing compo- 
nent is indeed critical for deriving column densities. From 
the upper panel of Fig. 3, allowing, conservatively, for 
3cr variations (about 10 for Ax 2 ), we find a huge range 
of values 14.6 < log 10 iV(HD) < 16.9, the lower end be- 
ing even larger than most of the Copernicus measure- 
ments. From the present data alone, the saturation of 
the lines prevents a more precise determination. However, 
one sees that the best estimates (Ax 2 ~ 0) are found for 
log 10 ./V(HD) around 16.1 (solid line in Fig. 2) which cor- 
responds, according to the bottom panel of Fig. 3, to an 
effective 6- value of about 1 km s _1 . 

3. DISCUSSION. 

Previous observations of molecules toward HD 73882 in- 
clude CH, CH+, CN, C 2 and CO (Gredel, van Dishoeck 
& Black 1993, 1994). Although the molecular fraction 2 
N(R 2 )/[N(H) + 2A(H 2 )] is quite close to that found to- 
ward C Oph (Snow et al. 2000), the CO/H 2 and CN/H 2 ra- 
tios, being very sensitive to the radiation field, differ 
greatly. These ratios are over 5 times larger for HD 73882 
and closer to values found for other dark clouds such as 
TMC1. The extinction curve for HD 73882 is peculiar 
compared to "normal" diffuse cloud curves, with a steep 
far-UV rise and a weak bump at 2200 A. Moreover, the 
value of R v =A v /i?B-v is 3.39, larger than the mean galac- 
tic value ~ 3 (Cardelli et al. 1989). 

Ground based mm- wave 12 CO reveal three components 
in emission separated by 3 and 2.3 km s _1 (i.e. 5.3 km 
s _1 for the extreme ones, nearly as for 13 CO; Gredel, van 
Dishoeck & Black 1994), not resolved with FUSE spectra if 
indeed they are really intercepted by the HD 73882 line of 
sight. Snow et al. (2000) show that Nal has many compo- 
nents. However, all diatomic molecules seen in absorption 
show only one component, including the C 2 observations 
recorded at high spectral resolution (R=10 5 ). These com- 
ponents agree in velocity. Therefore, several components 
might be present, but the possibility of a single molecu- 
lar absorbing region on the HD 73882 line of sight is not 
excluded. 

The present observation of this relatively high extinction 
target opens the possibility that, indeed, we are probing a 
cloud which could have reached the transition point where 
the reservoir of deuterium is primarily the HD molecule. 
We have modeled the photodissociation region with a den- 
sity of 350 cm~ 3 , as determined from the analysis of C 2 
excitation (Gredel, van Dishoeck & Black 1993), a stan- 
dard ultraviolet interstellar radiation field and a cosmic 
ionization rate £ = 5 x 10 s , considered as an aver- 
age (Abgrall et al. 1992; Le Bourlot et al. 1993; Roueff 
& Node-Langlois 1998). The atomic-to-molecular transi- 
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tion is assumed to be represented by a semi-infinite, plane 
parallel slab with the radiation field impinging on one 
side only of the cloud. The mechanisms involved in the 
photodissociation of HD are very similar to those for H2 , 
but self-shielding is much more efficient for the abundant 
H2 than for HD. Similarly to H2, HD is formed on grains. 
However, HD is also predominantly formed through the re- 
action H 2 + D + — > HD + H + , whenever substantial H 2 is 
present. 

The main results of this model are displayed in Fig. 4 
which shows the evolution of the fractional abundances 
of atomic hydrogen and deuterium, and their associated 
molecular species, as a function of the extinction A v . It 
is clear that the H to H2 transition occurs at significantly 
lower extinction than the D to HD transition (note that 
the C to CO transition would take place at even higher ex- 
tinction than HD). This reflects the smaller contribution of 
self-shielding in the photodissociation mechanism for HD. 
According to the extinction value for our target star, if 
there arc not a large number of velocity components, we 
might be dealing with an absorbing region where HD is the 
primary reservoir of deuterium. It also shows that even if 
several components are detected in atomic lines, they may 
be present in the H 2 lines, but not necessarily in the HD 
ones which are tracing a deeper region in the cloud. 

4. CONCLUSION. 

We have reported FUSE observations in the direction of 
the reddened star HD 73882. Several lines from the HD 



molecule are detected, together with many from H 2 . Un- 
der the assumption of a single absorbing molecular com- 
ponent, our best estimate for the HD/H 2 column density 
ratio from a simultaneous profile fitting of all detected lines 
is of the order of 10~ 5 . However, strictly speaking this is 
only an upper limit since the fit is very sensitive to the 
intrinsic width of the lines, most of them being saturated. 

This preliminary evaluation toward such a reddened star 
opens the possibility of probing regions sufficiently deep 
into interstellar molecular clouds for HD to be the main 
reservoir of deuterium atoms. In such regions, the mea- 
surement of the HD/H 2 ratio would thus become a new 
reliable and powerful tool for determining the interstellar 
D/H ratio in dense regions for which the simple relation 
(D/H) = 0.5 x HD/H 2 will apply. Recall that the possi- 
bility of multiple clouds for the present line of sight pre- 
vents claiming it is indeed the case. One may hope that 
a number of similar future FUSE observations will ulti- 
mately succeed, providing additional information on the 
structure of the observed sight-lines is available. 
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Fig. 1. — FUSE spectrum of HD 73882 from 1030 to 1120 A. Many atomic and molecular absorption lines are detected. Most of them are 
from molecular hydrogen (tick marks indicate those from J = 1 to J = 4). Lines from H2(J = 0) at 1037, 1049, 1063, 1077, 1092 and 1108 
A and from H2(J = 1) are strongly saturated. A CO band at 1088 A is also indicated. The gap around 1085 A is due to the gap in between 
the two micro-channel plates mounted in each LiF detectors. 

Fig. 2.— Plots of four B 1 S+-X 1 E+ electronic transitions of HD detected in absorption toward HD 73882 at 1042.85, 1054.29, 1066.27 and 
1105.91 A, corresponding respectively to J=5— >0, J=4— >0, J=3— >0 and J=0— >0 (the J=2— >0 and J=l— >0 lines are hidden by strong H2 lines). 
The actual FUSE wavelength absolute calibration being uncertain by few km s _1 , the wavelength scale has been shifted and the HD lines 
aligned at the same radial velocity arbitrarily set to km s . The solid line is the best fit obtained simultaneously for all detected HD lines 
and other identified absorptions. To underline the detection independently proven by the x 2 variation (see text and Fig. ??), the bottom 
panel shows the composite solution covering the first three HD lines, thus enhancing the S/N ratio (the line at 1105.91 A is blended with Cl 
and has been excluded in the summation). 

Fig. 3. — The upper panel is a plot of the A\ 2 of the spectral fits as a function of the HD column density. The lower panel gives the best 
estimates of ./V(HD) as a function of the effective b- value. Small b- values are slightly favoured by the A\ 2 . 

Fig. 4. — Comparison of the relative evolution of H, D, H2 and HD fractional abundances (relative to the proton density) with A v . We 
assume that a semi infinite plane-parallel cloud is exposed to an isotropic "standard" interstellar ultraviolet radiation field (scaling factor 
X=l). The cloud has a constant density riH=350 cm -3 , deduced from C2 observations toward HD 73882. £ is the cosmic ionization rate. For 
HD 73882, A v =2.44. The small changes in the HD and D curves beyond A v ~5 are due to slight variations in temperature within the cloud 
which in turn induce slight chemical variations. 
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